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Abstract On 15 January 2022, Hunga Volcano in Tonga produced the most violent eruption in the modern 
satellite era, sending a water-rich plume at least 58 km high. Using a combination of satellite- and ground-based 
sensors, we investigate the astonishing rate of volcanic lightning (>2,600 flashes min~') and what it reveals 
about the dynamics of the submarine eruption. In map view, lightning locations form radially expanding rings. 
We show that the initial lightning ring is co-located with an internal gravity wave traveling >80 m s~! in the 
stratospheric umbrella cloud. Buoyant oscillations of the plume's overshooting top generated the gravity waves, 
which enhanced turbulent particle interactions and triggered high-current electrical discharges at unusually 
high altitudes. Our analysis attributes the intense lightning activity to an exceptional mass eruption rate 

(>5 x 10° kg s~'), rapidly expanding umbrella cloud, and entrainment of abundant seawater vaporized from 
magma-water interaction at the submarine vent. 


Plain Language Summary The eruption of Tonga's underwater Hunga Volcano culminated on 

15 January 2022 with a giant volcanic plume that rose out of the ocean and into the mesosphere. This plume 
created record-breaking amounts of volcanic lightning observed both from space and by radio antennas on the 
ground thousands of kilometers away. We show that the eruption created more lightning than any storm yet 
documented on Earth, including supercells and tropical cyclones. The volcanic plume rose to its maximum 
height and expanded outward as an umbrella cloud, creating fast-moving concentric ripples known as gravity 
waves, analogous to a rock dropped in a pond. Point locations of lightning flashes also expanded outward in a 
pattern of donut-shaped rings, following the movement of these ripples. Optically bright lightning was detected 
at unusually high altitudes, in regions of the volcanic cloud 20-30 km above sea level. Our findings show that a 
sufficiently powerful volcanic plume can create its own weather system, sustaining the conditions for electrical 
activity at heights and rates not previously observed. Overall, remote detection of lightning contributed to 

a detailed timeline of this historic eruption and, more broadly, provides a valuable tool for monitoring and 
nowcasting hazards of explosive volcanism worldwide. 


1. Introduction 


An explosive eruption began on 19 December 2021 (coordinated universal time, UTC) near the remote islands 
of Hunga Tonga and Hunga Ha'apai in the South Pacific Ocean (Carn et al., 2022; Gupta et al., 2022). These two 
islands are small peaks on the rim of a much larger submarine caldera volcano known as Hunga Volcano (Taan- 
iela Kula, Tonga Geological Services, written communication, 10 September 2022). The explosive activity inten- 
sified on 13 January 2022 (UTC), followed by the climactic eruption on 15 January that sent a water-rich volcanic 
plume into the mesosphere (Carr et al., 2022; Podglajen et al., 2022; Proud et al., 2022). This event continues to 
push the boundaries of our understanding of how explosive volcanism impacts the broader Earth system. In addi- 
tion to significant local impacts on the Kingdom of Tonga (Purkis et al., 2023), the eruption created global-scale 
acoustic waves (Matoza et al., 2022), tsunamis (Omira et al., 2022), ionospheric and geomagnetic disturbances 
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(Astafyeva et al., 2022; Gavrilov et al., 2022), and warmed the climate due to water vapor injection (Sellitto 
et al., 2022). Additionally, Briggs et al. (2022) noted the high rates of volcanic lightning during this eruption, 
leading to key questions about the nature of electrification during such an extreme event. Here, we examine the 
evolution of the volcanic plume and its extraordinary lightning activity. Although we focus primarily on 15 Janu- 
ary, we also provide a data set of lightning characteristics from 13 to 14 January for future research. 


2. Data and Methodology 
2.1. Volcanic Plume Observations 


The volcanic plume on 15 January 2022 was characterized using visible and infrared observations from two 
geostationary satellites: GOES-17 (operated by the U.S. National Oceanic and Atmospheric Administration, 
NOAA) and Himawari-8 (Japan Meteorological Agency). All times are reported in UTC (local Tonga Time 
is UTC +13 hr). Our analysis accounts for the ~7 min difference between the scan start and the timing of the 
scan line over Hunga Volcano for both satellites. This time offset means that the 04:00 image does not record 
observations near the volcano until 04:07, which has important implications for establishing the chronology of 
the eruption. 


We retrieved stereoscopic cloud heights from the combination of nearly synchronous satellite images at 10-min 
intervals using 10.3-um band infrared brightness temperatures from GOES-17 and Himawari-8, as described in 
Podglajen et al. (2022). These height retrievals are accurate within +0.5 km and eliminate the effects of parallax 
related to the viewing angle of individual satellites. However, the analysis was limited to clouds >10 km above 
sea level (asl) and within a ~1,920 x 1,080 km grid around the volcano. Maximum plume heights were extracted 
within 100 km of Hunga Volcano for the first 7 hours of eruption and within 50 km thereafter. The lower-level 
plumes occurring early on 15 January were estimated using shadow lengths or minimum brightness temper- 
atures matched to atmospheric temperatures (Oppenheimer, 1998) with uncertainties of +3 km (Figure S2 in 
Supporting Information S1). Volcanic cloud areas were determined from the stereopairs from 04:17—07:07 and 
NOAA's GOES-17 based volcanic cloud detection and tracking product (Pavolonis et al., 2018) from 07:17—16:07 
(Van Eaton et al., 2023). We use the term umbrella cloud to refer to a density-driven, horizontal intrusion into the 
atmosphere (Sparks et al., 1997: Ch. 11). 


2.2. Lightning Detection and Analysis 


Several long-range systems detected the lightning from this eruption (Text S1, Table S1 in Supporting Infor- 
mation S1). First, we obtained data from the satellite-based Geostationary Lightning Mapper (GLM) onboard 
GOES-17. GLM uses pixel mapping to obtain the timing, location, flash area, and optical brightness of lightning. 
Next, we created a combined data set from three ground-based networks of radio antennas: (a) Vaisala's Global 
Lightning Data set (GLD360), and (b) the Earth Networks Total Lightning Network (ENTLN), which incorpo- 
rates data from (c) the World Wide Lightning Location Network (WWLLN). These systems detect radio waves 
emitted by impulsive discharges of electrical current, known as pulses, which occur intermittently throughout a 
lightning flash. Data include the timing, latitude, longitude, and estimated peak current of each lightning pulse. 
A lightning flash may consist of multiple pulses, which can be grouped together using temporal and spatial 
constraints (e.g., Cummins et al., 1998). For this study, pulses from the ground-based networks were combined 
into a single data set. Duplicates were removed if they occurred within 100 us and 5 km of a prior pulse. Then the 
pulses were clustered into lightning flashes serially by requiring the time between successive pulses to be no more 
than 100 ms and the total flash duration to be no more than 1 s. To be included in the flash, a pulse was required 
to be <10 km from at least one other pulse in the flash (see Text S2 and Figure S2 in Supporting Information S1). 


To minimize lightning contamination from non-volcanic storms, we identified nearby convective cells that did 
not originate from Hunga Volcano using GOES-17 multispectral imagery through NASA's Worldview platform. 
This process was relatively straightforward for 15 January, leading to the following constraints on volcanic light- 
ning: (a) within 250 km of the volcano for optical GLM observations, (b) within 300 km for ground-based obser- 
vations, and (c) constrained to 00:00-15:12 for all data. Although this paper focuses on the 15 January activity, 
we also include lightning data from 13 to 14 January 2022 to provide a complete data set for future researchers. 
We caution that there was substantially more complex storm activity on 13-14 January that merged with the 
volcanic clouds. We have not attempted to quantify or remove the non-volcanic lightning on those days, which 
will require careful attention in future work. 
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Figure 1. Satellite and volcanic lightning chronology of the 15 January 2022 eruption of Hunga Volcano in Tonga. Four 
phases of eruptive activity can be distinguished using umbrella cloud growth, maximum plume heights, and lightning rates. 
Lightning detected by ground-based networks is shown on the left axis (log scale) in terms of flashes min7! occurring within 
300 km (black line) and 20 km of the volcano (blue line). Optically detected flash rates (yellow bars) are shown on the right 
axis (linear scale) and include lightning detected within 250 km of the volcano. Timing of lightning rings indicated by two 
horizontal lines. Details of each eruption phase on 15 January are provided in Table S2 in Supporting Information S1. Heights 
of the tropopause (~17 km) and stratopause (~43 km) are shown as dashed lines on the plume height axis, given in km above 
sea level (asl). See Figure S3 in Supporting Information S1 for meteorological sounding data. 


Key issues related to the configuration of each detection system are described in the Text S1, Figure S1, and Table S1 
in Supporting Information S1. One limitation is that unusually high flash rates saturated the GLD360 sensors from 
04:30 to 06:00 on 15 January. During this time the system prioritized detection of higher-current lightning, which 
meant that many of the low-current pulses were lost. These effects show up clearly in the lightning locations within 
20 km of the volcano as an apparent dip in proximal flash rates (Figure 1, Figure S1 in Supporting Information S1). 


3. Results and Analysis 
3.1. Chronology of Eruptive Activity on 15 January 


The timelines of plume heights, umbrella cloud growth, and lightning flash rates reveal distinct phases of eruptive 
activity on 15 January (Figure 1, Table S2 in Supporting Information S 1). Significant airborne emissions lasted for at 
least 11 hr. Phase 1 began with an initial volcanic plume visible in satellite images by 02:57. The plume rose briefly to 
~15 km (based on shadow lengths) before returning to lower altitudes by 03:20 and persisting through 03:57. There 
was no lightning detected yet. A sudden increase in the plume's width occurred in the 04:07 satellite images. Vergoz 
et al. (2022) identified the onset of long-range infrasound at 04:00, and we suggest that a step-change in the eruption 
dynamics may have led to a larger, higher velocity plume at this time. Phase 2 encompasses the most intense stage 
of the eruption from 04:00-07:11, with mass eruption rates (flux of magma, magmatic gas, and vaporized seawater) 
exceeding 5 x 10° kg s~! from 04:15-04:57 based on the rate of umbrella growth (Text S3, Figure S4 in Supporting 
Information S1). The plume rose to a maximum overshoot of 58 km and expanded into at least two atmospheric 
levels, as noted by Carr et al. (2022) and Gupta et al. (2022). The upper umbrella intruded into the stratosphere 
(reaching 30-40 km and higher) while a second, lower cloud spread out at the level of the tropopause (but with cloud 
tops reaching ~20 km). A renewed phase of eruptive activity began at 07:46 with the return of proximal lightning 
and an overshooting top 39 km high by 08:47 (Figure 1). Phase 4 represents lower-level activity with an overshoot 
of 28 km at 13:37. Based on the end of proximal volcanic lightning, we estimate the eruption ended around 15:12. 


3.2. Comparison of Lightning From Radio and Optical Systems 


One striking observation is the gap in optically detected lightning just as plume heights and lightning rates from 
ground-based radio frequency systems reached their peak (Figure 1, Figure S5 in Supporting Information S1). 
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The earliest lightning on 15 January was detected by ground-based systems at ~04:09 as the plume ascended 
rapidly from 5 to 23 km asl. GLM observed 22 optical flashes from 04:14-04:25, which became gradually 
smaller and dimmer as the plume continued to rise (flash size decreased from 8,243 to 4,487 km? and optical 
energy decreased from 1,543 to 86 fJ). The lightning became undetectable by GLM as the plume rose above 
~30 km asl despite the dramatic increase in flash rates from ground-based networks (Figure 1; eruption Phase 
2). These observations suggest two different possibilities. One is that lightning may have occurred only below 
~30 km and the optical energy could not penetrate the overlying cloud (Schultz et al., 2020). Alternatively, it is 
possible that lightning did continue throughout the entire cloud, but at heights above ~30 km the flashes were 
too small or optically dim to be detected by GLM (cf. Cummins, 2021; Rison et al., 2016). However, we do not 
favor this explanation because the very high peak currents detected from 04:30-05:40 (Figure S1 in Supporting 
Information S1) suggest that any traditional lightning occurring near the cloud top would have been optically 
bright (Quick & Krider, 2013). We explore these concepts further in Section 4.3. 


By 05:37, the eruptive intensity weakened enough that stratospheric winds began to shear off the upper umbrella 
and reveal parts of the lower cloud below (see Figure 2 and Figure S3 in Supporting Information S1). Yet optical 
detections did not resume until 4 hours later (at 08:32) during a phase of renewed activity (phase 3a; Figure 1). 
Lightning could only be optically detected in regions of the cloud below ~30 km during this phase, determined 
by mapping the GLM flashes onto stereoscopic cloud height retrievals (Figure 2, Figure S5 in Supporting Infor- 
mation S1, Movie S1). The greatest number of optically detected flashes occurred from 12:17—15:21 once the 
upper umbrella drifted away from Hunga Volcano and exposed the core of the volcanic plume directly above the 
vent. The size of GLM flashes also reached their peak during this time (~20,000 km?) even though cloud tops 
were much lower than earlier in the eruption (17—28 km asl; Figure 2). 


3.3. Gravity Waves and Lightning Rings in the Expanding Umbrella Cloud 


It is common for strong thunderstorms, wildfire plumes, and volcanic plumes to produce an overshooting top 
that oscillates around the level of neutral buoyancy, sending out concentric gravity waves like ripples from a 
rock dropped into a pond (e.g., Apke et al., 2018; Bonadonna & Phillips, 2003; Holasek et al., 1996; Luderer 
et al., 2007). The Hunga eruption took this process to the extreme. Its overshooting dome was more than 100 km 
across and rose to 58 km asl by 04:37, then dropped to 30 km asl over the next 10 min (Figure 1, Figure $4 in 
Supporting Information S1). By 04:47, the resulting gravity wave was >10 km high from peak to trough and 
moved outward at speeds >80 m s~! (Figure S6 in Supporting Information S1). The vertical motion of this 
gravity wave superimposed on the lateral expansion of the density-driven umbrella cloud, which moved outward 
at rates of 60-90 m s~! from 04:27 to 04:47. For context, winds this strong at ground level (>130 miles hr!) 
would uproot large trees and cause catastrophic damage. In map view, the ground-based lightning locations reveal 
successive, radially expanding rings of electrical activity. The overlay of stereoscopic cloud heights and lightning 
locations show that the first and largest ring is co-located with the leading edge of the first and largest gravity 
wave in the upper umbrella cloud (see 04:16—04:52 in Figure 2, Movie S1). 


4. Discussion 


4.1. What Caused the Lightning Rings? 


The initial lightning ring concentrated at the leading edge of an enormous gravity wave in the stratospheric 
umbrella cloud (Figure 2, Movie S1). This feature-matching suggests the lightning ring occurred predominantly 
in the upper umbrella rather than in the lower-level portion of the cloud that spread out near the tropopause. Light- 
ning rings (also known as lightning holes) have been observed in severe storms before (Calhoun et al., 2013, 2014; 
Payne et al., 2010), but never associated with convective gravity waves, nor on this scale (~140 km radius for 
Hunga's largest ring compared to ~10 km in severe storms). It is possible that a distinct mechanism governed 
ring formation in this case. 


We propose that the lightning radiated outward in a concentric pattern because it was modulated by turbulence and 
air pressure changes induced by gravity wave propagation. Updrafts associated with the transverse (i.e., vertical) 
motion of the gravity waves—combined with horizontal expansion of the density-driven umbrella region—would 
have provided a source of particle motion to organize and separate regions of charge. Numerical simulations 
have shown that gravity waves from overshooting tops in wildfire plumes can induce turbulence and sharp gradi- 
ents in air pressure and particle concentration (Luderer et al., 2007). Hunga's fast-moving gravity waves likely 
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Figure 2. Maps of volcanic plume and lightning development on 15 January 2022, with times shown in UTC. Grayscale gives stereoscopic cloud heights, blue dots 
show lightning flashes detected by ground-based radio frequency networks over the following minute, and purple-yellow color scale shows optically detected lightning 
from the GLM sensor. (*) indicates frames with optically detected lightning. At least four distinct lightning rings occur from 04:16 to 05:51 followed by a final ring 
from 08:38-08:48. The initial and most prominent ring (visible in the first four frames) concentrated at the leading edge of a gravity wave within the upper umbrella 
cloud. Pink circles outline the lightning ring in two frames, showing an (average) expansion rate exceeding 60 m s~!. Westward advection of the upper umbrella starts 
to reveal a lower level cloud by 05:37. White dashed polygons outline the lightning locations, showing their westward movement with the stratospheric umbrella cloud. 


Local islands are outlined in black. 
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segregated the particles by size (and thus, electric charge) depending on how they coupled to the flow. These 
processes seem to have preferentially occurred at the wave front, creating lightning in a circular pattern. As the 
wave moved through, it triggered enough lightning to neutralize the charge locally, explaining the lightning 
“hole” from 04:37 to 04:47 (Figure 2). Yet 2 minutes later, the gap begins filling in, apparently from the ring 
inward (Movies S1, S2)—we do not currently have a satisfactory explanation for this phenomenon. Future studies 
of the mechanisms behind these unique observations would benefit from a three-dimensional, non-hydrostatic 
model of volcanic plume electrification. However, such a model does not yet exist (cf. Textor et al., 2006). 


4.2. Why Was the Lightning Activity So Intense? 


Our analysis demonstrates that the 15 January eruption produced 2,615 flashes min™! at its peak intensity (04:53), 
representing the most intense electrical storm ever detected by global networks (Table S3 in Supporting Informa- 
tion S1). This peak lightning rate is significantly higher than the next most lightning-rich case study in our analysis 
(cf. 993 flashes min~!, Table S3 in Supporting Information S1). What made this eruption so intensely electrified? 
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Figure 3. Conceptual sketch of dynamic and microphysical processes in the volcanic plume on 15 January 2022. Plume 

top morphologies are based on transects through the stereoscopic cloud height retrievals (Figure S6 in Supporting 
Information S1), whereas the underlying structure is inferred. Blue shading indicates the mixed phase zone, where liquid 
water, cloud ice, graupel, and ash coexist—this is where we expect most of the lightning to occur. (a) Represents the 
developing plume before lightning was detected (before 04:09 UTC); (b) represents the plume at its maximum height >50 km 
when oscillation of the overshooting top sent out gravity waves through the upper umbrella (04:37-05:27); (c) represents the 
waning eruption as the upper umbrella detached and drifted west (after 05:37). 


Previous studies have shown that volcanic plumes become electrified in two ways: (a) silicate charging from ash 
particle breakup and collisions and (b) ice charging, which is analogous to regular thunderstorm electrification 
(Mather & Harrison, 2006; Cimarelli et al., 2022, and references therein). The submarine location of Hunga's 
vent allowed the eruption to vaporize vast amounts of seawater, transporting an unprecedented 50-150 Tg of 
water vapor into the stratosphere (Millan et al., 2022; Vomel et al., 2022). In general, “wet” eruption plumes start 
at magmatic temperatures (~1000°C) and cool during explosive mixing with surrounding water (Koyaguchi & 
Woods, 1996). The high initial velocity and heat budget of major volcanic eruptions allow them to penetrate much 
higher into Earth's atmosphere compared to severe thunderstorms (Herzog & Graf, 2010; Sparks et al., 1997; 
Woods, 1995). Ice charging occurs when a plume creates strong updrafts through a mixed phase region containing 
liquid and frozen hydrometeors. These conditions are much more likely to arise in large volcanic plumes created 
by magma-water interaction (e.g., Behnke et al., 2013; Prata et al., 2020; Van Eaton et al., 2012, 2020, 2022) 
because the added water leads to high concentrations of liquid and frozen hydrometeors in the upper plume. In 
fact, large-eddy simulations show that it is theoretically possible to sustain liquid water in the warm interior of an 
eruption column well into the stratosphere (Van Eaton et al., 2012; their Figure 3). Therefore, our observations 
of optically detected lightning up to 20-30 km asl (Section 3.2, Figure S5 in Supporting Information S1) suggest 
that the lightning-producing, mixed-phase region of Hunga's plume may have also extended up to ~30 km high. 
Figure 3 shows a conceptual sketch of a mixed-phase microphysical zone in the warmer core of the plume, carry- 
ing liquid water, icy hydrometeors, and ash into the mid- to upper stratosphere. 


However, it is not just the plume's water content that produced the extraordinary storm. Volcanic lightning rates 
scale with eruptive intensity—taller plumes and broader umbrella clouds produce more flashes (e.g., Behnke 
et al., 2013; Smith et al., 2021; Van Eaton et al., 2022). The Hunga eruption shows positive relationships among 
flash rate, plume height, and umbrella expansion through time on 15 January (Figure S4 in Supporting Informa- 
tion S1). In supercell thunderstorms, flash rates are correlated with both the updraft mass flux and the flux of 
precipitating ice (graupel) due to increased hydrometeor collisions (Deierling et al., 2008; Kulman et al., 2006). 
Turbulence creates a clumpy distribution of charge that supports higher flash rates due to the increase in high 
electric field regions at the boundaries of charge pockets (Bruning & MacGorman, 2013). Flash rates have been 
shown to correlate with turbulent kinetic energy (Behnke & Bruning, 2015; Bruning & MacGorman, 2013). 
Therefore, it is likely the high mass flux in Hunga's volcanic plume (>5 x 10? kg s~! by our calculations) encour- 
aged more turbulent inter-particle collisions that enhanced charge generation. 


4.3. Evidence of Stratospheric Lightning Initiation? 


Optically bright lightning was detected by the GLM sensor in regions of the volcanic plume as high as 20-30 km 
asl (Movies S1, S2, Figure S5 in Supporting Information S1). Although there are large uncertainties on these 
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estimates, the possibility of lightning at such high altitudes is intriguing. Previous studies have only detected 
lightning initiation processes up to ~18 km asl in overshooting cloud tops (DiGangi et al., 2016; MacGorman 
et al., 2017). In Earth's atmosphere, it becomes easier for electrical breakdown to occur at higher altitudes, but 
more difficult to create the hot plasma channels (leaders) that characterize lightning (Dwyer & Uman, 2014). 
For example, electrical phenomena known as blue jets and gigantic jets begin as lightning leaders propagating 
out of a thundercloud top, but only the cold plasma (streamer zone) propagates into the mid stratosphere and 
beyond (Krehbiel et al., 2008). The theoretical limit on this streamer-to-leader transition in Earth's atmosphere 
is 20-23 km asl, due in part to the very low air densities at high altitude, which inhibit the heating of air (Boggs 
et al., 2022; van der Veld et al., 2019). The transition is dictated by the current and radius of the leader, among 
other characteristics (da Silva and Pasko, 2013). Therefore, our evidence for optically bright flashes (i.e., leader 
forming activity) detected by GLM up to ~30 km asl would require an unconventional process. One possible 
mechanism is that strong pressure anomalies created by the fast-rising plume (e.g., Oberhuber et al., 1998, their 
Figures 2—3) and its overshooting top (e.g., Luderer et al., 2007, their Figure 8) created locally higher air pressures 
that supported conditions for leader formation at extreme heights. 


In addition to high-altitude lightning, the eruption may have produced a special class of electrical discharges 
known as narrow bipolar events (NBEs). NBEs are distinct from traditional lightning because they consist 
entirely of cold plasma (streamers). NBEs are believed to be responsible for initiating lightning, though they 
often occur in isolation (Rison et al., 2016). They have historically been associated with severe convection and 
overshooting tops in thunderclouds (Karunarathna et al., 2015; Wu et al., 2012, 2013) and produce mainly blue 
light (~337 nm) in the visible regime. Blue light is not detectable by the satellite-based GLM sensor (~777 nm) 
and could explain why there were no optical detections during the peak of the eruption—even if NBEs occurred 
right at the plume top, they would be invisible to GLM. Likewise, the exceptionally high peak currents detected 
by the ground-based networks during lightning ring activity (from 4:30-5:40, Figure S1 in Supporting Informa- 
tion S1) are consistent with NBEs (Le Vine, 1980; Smith et al., 1999). Some of these high-current discharges 
appear to follow the stratospheric cloud as it drifts downwind (white dashed outlines, Figure 2). Furthermore, our 
observation that traditional flash clustering algorithms do not seem to apply to these lightning data (Text S2 in 
Supporting Information S1) provides support for isolated electrical discharges like NBEs. Isolated, blue stream- 
ers have been photographed in a volcanic plume before (Van Eaton et al., 2022). Therefore, we remain open to the 
possibility that some (or all) of the “lightning rings” could be attributed to NBEs rather than lightning in the strict 
sense. This scenario would still imply an upper limit of ~30 km for traditional (optically bright) lightning in the 
plume, as detected by GLM. And since the highest reported altitudes of NBEs are ~20 km (Wu et al., 2012), the 
overall evidence for stratospheric electrical activity in Hunga's plume suggests that volcanic eruptions may create 
conditions for lightning initiation (i.e., leader and/or NBE formation) at unprecedented altitudes. 


5. Conclusions 


Our volcanic lightning and plume analysis reveals at least four phases of eruptive activity from 02:57 to 15:12 
on 15 January 2022. In particular, the rates of proximal lightning (within 20 km of Hunga Volcano) resolve 
the timing and intensity of volcanic plume formation, augmenting satellite observations when the region was 
obscured by the upper-level cloud for several hours (especially ~10:00—15:00, Figures 1 and 2). These results 
frame an eruptive timeline for this globally important, yet remote, eruption, with implications for monitoring 
volcanic hazards using real-time lightning data in the future. We expect the eruption chronology to evolve as more 
details from near-source observations, geological deposits, and eyewitness reports come to light. 


This eruption produced the largest lightning rings ever observed (up to ~140 km radius). We show that the 
initial lightning ring matches the movement of a large-scale gravity wave in the upper umbrella cloud, formed 
by dramatic oscillation of the overshooting top (Figure 2, Figure S6 in Supporting Information S1). Results also 
indicate that the 15 January eruption produced the highest lightning rates yet documented (>2,600 flashes min~'). 
We explain the intense electrification as a function of the exceptionally high eruption rate (at least 5 x 10° kg s7! 
including magma and external water), fast-growing umbrella cloud, and incorporation of seawater from the 
submarine vent. Comparison of different lightning detection technologies from spaceborne optical measurements 
and ground-based radio antennas constrain the height of lightning-producing regions of the volcanic plume. 
Based on the occurrence of optically bright flashes at high altitudes, we infer that the mixed-phase region of the 
volcanic plume extended up to ~30 km high, implying that liquid water injected deep into the stratosphere. Future 
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numerical simulations of the volcanic plume will be able to test the evolution and longevity of this mixed-phase 
microphysical region. Overall, our observations suggest that volcanic plumes can create the conditions for light- 
ning initiation well outside the range of meteorological thunderstorms previously observed in Earth's atmosphere. 
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GOES-17 and Himawari-8 data are accessible through https://registry.opendata.aws/noaa-goes/, https://regis- 
try.opendata.aws/noaa-himawari/, and NASA Worldview (https://worldview.earthdata.nasa.gov/). Stereoscopic 
cloud heights are downloadable from https://science-data.larc.nasa.gov/LaRC-SD-Publications/2022-11-22- 
001-KMB/. The data set of lightning and plume characteristics is downloadable from: https://doi.org/10.5281/ 
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